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Turbulent Oxygen Flames in Type la Supernovae 
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ABSTRACT 

o. 

In previous studies, we examined turbulence-flame interactions in carbon- 
burning thermonuclear flames in Type la supernovae. In this study, we consider 
turbulence-flame interactions in the trailing oxygen flames. The two aims of 
the paper are to examine the response of the inductive oxygen flame to intense 
levels of turbulence, and to explore the possibility of transition to detonation 
in the oxygen flame. Scaling arguments analogous to the carbon flames are 
presented and then compared against three-dimensional simulations for a range of 
'-pH i Damkohler numbers (Da^) at a fixed Karlovitz number. The simulations suggest 

that turbulence does not significantly affect the oxygen flame when Dai6 < 1, and 
the flame burns inductively some distance behind the carbon flame. However, for 
Dai6 > I? turbulence enhances heat transfer and drives the propagation of a flame 
that is narrower than the corresponding inductive flame would be. Furthermore, 
burning under these conditions appears to occur as part of a combined carbon- 
oxygen turbulent flame with complex compound structure. The simulations do 



> 



not appear to support the possibility of a transition to detonation in the oxygen 
(S| ■ flame, but do not preclude it either. 

Subject headings: supernovae: general - - white dwarfs - - hydrodynamics - 
nuclear reactions, nucleosynthesis, abundances - - conduction - - methods: nu- 
merical — turbulence — distributed flames 
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1. INTRODUCTION 

A major uncertainty in the modeling of Type la supernovae (SN la) is the phys- 
ical process whereby a subsonic deflagration tran sitions to a detonation. Such a tran 



sition seems to be required by the observations (jKozma et al.l 120051 ; iHoflich et al.l 11995 
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Mazzali et all 120071 ; iKasen et al.l 120091 ) , at least within the context of the popular "sin- 



gle degenerate" Chandrasekhar-mass mod el. Pr evious papers (e.g. iKhokhlov et al.l 119971 ; 



Niemeyer &: Woosleyl 119971 ; IWoosley et al.l 120091 ; I Asp den et al.l |2010| ) have focused on the 
possibility of a transition to detonation in carbon-rich material as the deflagration enters the 
"distributed burning regime" where burning becomes slow enough for turbul ence to disrupt 
the flame, mixing hot ash and cold fuel. Other papers (e.g. iPlewa et al.l 120041 ) have explored 
the possibility that detonation may be mechanically induced by the collision of burning waves 
near the surface of the white dwarf. The present situati on is inconclusive. T he collisions may 
not be strong enough to robustly cause a detonation (jRopke et al.l 120071 ) and the amount 
of turbulence require d for a spontaneous detonation in the distributed regime is quite large 
rtWooslevet"ai]l2009l ). 



Here we consider a third possibility - that the necessary carbon detonation actually be- 
gins as an oxygen det onation in a hybrid flame. This possibility has been recently considered 
( IWoosley et al.ll2010l ) in a one- dimensional study. Carbon burning produces oxygen-rich ash 
that still contains a large potential reservoir of nuclear energy. The oxygen ash is produced, 
for a given fuel density, at a constant temperature that gradually rises as a result of oxygen 
burning, until, finally, a silicon-rich composition is produced. As a result of turbulence, 
this oxygen layer, which we shall refer to as an oxygen "flame" , is broadened and islands of 
nearly isothermal conditions are produced. Here we flesh out those one-dimensional results 
in a series of three-dimensional simulations. 



In two previous three-dimensional studies, I Asp den et al.l (j2008al . |2010| ) (henceforth Pa- 
pers I and II), turbulence-flame interactions in carbon-burning flames were examined at 
small and large scales, respectively. In Paper I, it was shown that once the turbulence was 
sufficiently strong, the mixing of fuel and heat was driven by turbulent mixing instead of 
thermal diffusion. This resulted in a categorically different kind of flame, which was referred 
to as a distributed flame. Paper II extended these small-scale st udies to (more rea listic) 
larger length scales, where scaling relations based on the theory of iDamkohlerl (119401 ) were 
predicted to reach a limiting behavior, resulting in a so-called "A-flame". 

In this paper, we focus on turbulence-flame interactions in the trailing oxygen flame, 
which are expected to be significantly different than in the carbon flame, due to the inductive 
nature of the oxygen burning. The specific question we address is whether turbulence can lead 
to a greatly extended oxygen-rich region that might have properties suitable for detonation. 
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2. THEORETICAL DESCRIPTION 

We first recap the scaling relations for turbulent carbon flames from Paper II, and 
then present theory describing possible modes of burning in turbulent oxygen flames. We 
fix the fuel to be carbon-oxygen at a particular density and temperature, and then choose 
turbulent conditions based on the flame properties of the fuel. Since the oxygen burning 
time scale is determined by the temperature resulting from carbon burning at the given 
density, the width of the oxygen flame is most sensitive to that density and the turbulence 
properties. For typical turbu lent conditions i n the supernova and an initial composition of 



40% carbon and 60% oxygen, IWoosley et al.l (120 lOl ) find that the density of greatest interest 



is P12 = 2.5 x 10 7 g cm -3 (we use p and T for total density and temperature the suffix 12 
denotes conditions before the carbon burns). This gives a post-carbon-flame temperature and 
density of approximately T\q = 3.14 x 10 9 K and p±e = 1.69 x 10 7 g/cm 3 , respectively (here 
suffix 16 denotes conditions after carbon burning but before oxygen has burned). Under 
these conditions, the oxygen flame ha s an inductive burning time scale of approximately 



n 6 = 0.016 s, see IWooslev et al.l (120 lpf ) 



Having fixed the fuel conditions, the two parameters that can be varied are the rms 
turbulent velocity fluctuation u and the integral length scale /. Turbulent premixed flames 
are characterized through Karlovitz and Damkohler numbers 

^ 2 U% S L l 

where sl and II are the laminar flame speed and width, respectively. These quantities 
represent the ratio of turbulent time scales at the Kolmogorov and integral length scales, 
respectively, and are two dimensionless quantities that represent the parameter space. As 
in Paper II, we focus on a fixed Ka^, corresponding to fixing the energy dissipation rate 
e* = u 3 /l of the turbulence in the star. For fixed Ka^ and e*, u = (e*/) 1 / 3 , it can be shown 
that D&l oc Z 2 / 3 . Therefore, the parameter space is one-dimensional and can be represented 
equivalently by either Da^ or /. 

We assume that the Karlovitz number is constant and sufficiently high to obtain a 
distributed carbon flame, see Paper I. The turbulent carbon flame properties will then depend 
on its turbulent nuclear burning time scale t^ 2 (note the superscript T differentiates the 
turbulent from laminar burning time scales, again see Paper I), and the properties of the 
turbulence, specifically the integral length scale / (recall that the turbulent intensity at 
fixed Karlovitz number is determined by u = (e*!) 1 ^ 3 ). Following Paper II and [Damkohler 



(jl940l ). by analogy with laminar flames, the turbulent flame speed sf 2 an d width l\ 2 can be 
expressed in terms of r^ and a turbulent diffusion coefficient T>t (not to be confused with 
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the Damkohler number Da) as 



s i2 = \I^Fi and l 12 = \/ t> tt^ (2) 



'12 



respectively. These relations only hold when the time scale of the turbulent eddies is shorter 
than the turbulent nuclear time scale of the carbon fuel, i.e. for Da 12 = t/t^ :Sl, where r is 
the turbulence time scale r = l/u. Taking a simple approximation T>t = ul, the turbulent 
flame speed and width were both shown to be proportional to Da 12 when Da 12 &1. For 
Da 12 £l, the turbulence can no longer broaden the flame, and the limiting A- flame behavior 
is reached (see Paper II), with local turbulent speed s A 2 and width A12 that depend on e* 
and t1 2 only, according to the relations 



A _ . L*^T „„a \ _ 4/ /*^r 3 



s i2 = \J e * r i2 and A u = Y£*rS , (3) 

respectively. Note that s A 2 and A12 are both constant. The turbulent flame speed and width 
can therefore be written as 

g i2 = l w = [ Da i2 for Daf 2 ^l, 
s i2 ^12 I 1 otherwise. 

We now apply this theoretical approach to oxygen flames, where there are three potential 
modes of burning: inductive, turbulent or A-flame. Each mode will have a corresponding local 
flame width and speed, which will be denoted (Zi6, s\q), (/f 6 , sf 6 ), and (Z A 6 , s A 6 ), respectively. 

The inductive mode is the simplest and is considered first. In a frame of reference where 
the carbon flame is stationary, the incoming fluid speed is equal to the flame speed Uq = s\i- 
The resulting oxygen flame has a width equal to /i6 = uqTi§, where the time taken for the 
oxygen to burn at a given density and temperature is t\q. In the presence of turbulence, 
the carbon flame speed is enhanced, but the oxygen flame remains slaved to ca rbon flame 



and li ft = MpT i fi, only wit h u = sf 2 . In the large-scale turbulence limit (see [Damkohler 



( 1l940l ); iPeterd ( 1l999l . 120001 )) . the turbulent carbon flame speed will be close to the turbulent 
intensity, and so we can take uq = (3ii, where /3 should be expected to be order unity, but can 
be as high as three, accounting for fluctuations and the density jump across the carbon flame. 
This defines si 6 = /3u and l w = (3utiq for a turbulent oxygen flame burning inductively. 

Defining Karlovitz and Damkohler numbers for oxygen flames 

Ka 16 = ^TT' and Dai 6 = ~1— > ( 5 ) 



6 16 t Ut 16 

reveals an interesting difference from carbon flames. Using s 16 = /3u, it can be shown that 
/3 2 Ka^ 6 Da 16 = 1, which means that the parameter space for oxygen flames is one-dimensional. 
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Therefore, under the assumption that the carbon flame is in the large-scale turbulence limit 
(so s 16 = /3u), the behavior of oxygen flames can be classified through the Damkohler number 
Dai6 alone. Note that for a fixed fuel, all of the relevant Damkohler numbers are constant 
multiples of each other, e.g. Da 16 = aDa[ 2 , where o = t^/ti6 is a constant. Also note that 



both si6 and lie can be shown to be proportional to Da 



1/2 

16 



If turbulent mixing can drive the flame, similar to the behavior in Paper II, scaling 
relations for turbulent flame speed and width can be predicted in terms of the oxygen burning 
time scale and the turbulent diffusion coefficient V T as s± 6 = ^/V T /ri e and lj 6 = y/VrTiQ. 
These scaling relations should only be expected to be possible for low values of the oxygen 
Damkohler number, i.e. Dai6 — ^1- O n the other hand, for Dai6^1, it may be possible to 
produce an oxygen A-flame, where the flame speed and width would be s^ 6 = y/e*TiQ and 



Ai6 = y/ £ * T w 3 i respectively As above, the turbulent oxygen flame speed and width can 
predicted to be 

s ie __ *i6 f Da ie for Daie^l, 

s i6 ^16 1 1 otherwise. 

The simulations in this study correspond to full-star conditions where it* = 2 x 10 7 cm/s 


be 

(6) 
on 



an integral leng th scale of L * = 1 x 10 6 cm, giving an energy dissipation rate of e* = 8 x 
10 15 cm 2 /s 3 , fsee lRopkdl2007[ ). Carbon fuel (40%) at p 12 = 2.5xl0 7 g/cm 3 andT 12 = 6xl0 8 K 



burns to pi$ = 1.69 x 10 7 g/cm 3 and T 16 = 3.14 x 10 9 K, which has an inductive time scale for 
oxygen of 0.016 s. This gives an oxygen A-flame speed and width of s^ 6 = 1.13 x 10 7 cm/s, 
and A 16 = 1.81 x 10 5 cm, respectively. 

Figure [1] depicts the scaling relations for the different turbulent flame widths as a func- 
tion of Damkohler number. Recall Dai6 c* ^ 2//3 > so this can be thought of as a function 
of integral length scale, shown by the thick black line. The red lines show the normalized 
inductive flame width he/Xie (solid for (3 = 1, and dashed for j3 = 2). The blue line shows 
the normalized turbulent width /f 6 /Ai6 = Dai6 for (Da^sl) and /f 6 /Ai6 = 1 for (Dai6<T), 
if turbulent mixing drives the flame. The red circles correspond to the simulations that will 
be considered in this study, specifically Dai 6 = 1/3, 1, and 3, with (3 = 1 and 2, and will be 
described in detail below. 

The black diamond highlights the A-flame thickness at Dai6 — 3, which is smaller than 
the corresponding inductive flame widths. This means that the possibility of an oxygen 
A-flame is particularly interesting as it would be an example of turbulence giving rise to a 
flame that is narrower than its zero turbulence counterpart. This is counterintuitive, as one 
expects turbulence to broaden interfaces. However, in this case, turbulence acts to enhance 
heat transfer, allowing the flame to burn more rapidly; turbulence can mix hot ash with cold 
fuel more rapidly than the fuel is heated through inductive burning. 
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SIMULATION DESCRIPTION 



As in Papers I and II, we use a low Mach n umber hydrodyn amics code, adapted to 
the study of thermonuclear flames, as described in lBell et al.l (120041 ) . The advantage of this 
method is that sound waves are filtered out analytically, so the time step is set by the bulk 
fluid velocity and not the sound speed. This is an enormous efficiency gain for low speed 
flames. We note that in all simulations presented here, the Mach number remains below 
0.1 (usually by an order of magnitude), and so c ompressibility effect s are c onsidered to be 
negligible. The reactions rate here are tak en fromlCaughlan fc Fowler! (119881 ) with screening. 
The conductivities are those reported in iTimmesI (J2000l ). a nd the equation of s tate is the 
Helmholtz free-energy based general stellar EOS described in lTimmes fc Swestyl (120001 ) . We 
note that we do not utilize the Coulomb corrections to the electron gas in the general EOS, 
as these are expected to be minor at the conditions considered. 

The non-oscillatory finite-volume scheme employed here permits the use of implicit 
large eddy simulation (iles). This technique captures the inviscid cascade of kinetic energy 
through the inertial range, while the numerical error acts in a way that emulates the dissipa- 
tive physical effects on the dynamics at the grid scale, without the expense of resolving the 



entire dissipation sub r ange. An overview of the technique can be found in iGrinstein et al. 



(120071 ). lAspden et al.l (j2008bl ) presented a detailed study of the technique using the present 
numerical scheme, including a characterization that allowed for an effective viscosity to be 
derived. Thermal diffusion plays a significant role in the flame dynamics, so it is explicitly 
included in the model. Species diffusion is significantly smaller, so it is not included explic- 
itly, but will be subject to numerical diffusion, which can be considered to ha ve an effective 



unity S chmidt number and exhibit the same behavior observed for viscosity in lAspden et al. 
(l2008bh . 



The turbulent velocit y field was maintained using the forcing term used in Papers I, II 
and lAspden et al.l (j2008bl ). Specifically, a forcing term was included in the momentum equa- 
tions consisting of a superposition of long wavelength Fourier modes with random amplitudes 
and phases. The forcing term is scaled by density so that the forcing is somewhat reduced 
in the ash. This approach provides a way to embed the flame in a turbulent background, 
mimicking the much larger inertial range that these flames would experience in a type la 
supernova, without the need to resolve the lar ge-scale convective motions that drive the tur- 
bulent energy cascade. lAspden et al.l (j2008bl ) demonstrated that the effective Kolmogorov 
length scale is approximately 0.28Ax, and the integral length scale is approximately a tenth 
of the domain width. 

Figure [2] shows the simulation setup. The simulations were initialized with oxygen 
fuel in the lower part of the domain and sulphur ash in the upper part, resulting in a 
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downward propagating flame. A high-aspect ratio domain was used, with periodic lateral 
boundary conditions, and outflow at the upper boundary. Due to the huge disparity in 
widths of carbon and oxygen flames, we are not able to capture both accurately in the 
same simulation. Therefore, we use carbon post-flame conditions and only burn oxygen. To 
recreate these condition appropriately, specifically the warm oxygen in a post carbon-flame 
mean velocity, we need to work in a moving frame of reference. Consequently, unlike papers I 
and II, a mean inflow was specified at the lower boundary to replicate the desired conditions 
so that an inductive oxygen flame could develop properly. Using such an inflow velocity in 
conjunction with the forcing term used to maintain the turbulent velocity field requires some 
care. It is possible to use a turbulent inflow velocity, but we have opted not to take such an 
approach. Instead we specify a uniform inflow and use just the forcing term itself to produce 
turbulence. We found that this gave satisfactory results, provided fi £2. 

The conditions of particular interest are when Dai6^1. Note from figure [1] that the 
integral length scale is expected to be much larger than all of the turbulent flame widths 
under these conditions. Given that the integral length scale is approximately one-tenth of 
the domain width, this means that the oxygen will burn extremely close to the inlet, and 
certainly before the turbulence has become well-developed. To account for this, the inflow 
velocity, temperature and density were synthetically altered (slower, cooler and more dense, 
respectively) to delay oxygen burning and ensure that the flame burned approximately half- 
way through the domain, while maintaining conditions close to the carbon post-flame. 

Oxygen Damkohler numbers of 1/3, 1 and 3 were simulated to capture the potential 
transition. The aim is to detect the mode in which the oxygen flame is burning, specifically, 
what is the local turbulent oxygen flame width. However, it is difficult to measure a local 
turbulent flame width directly, and the widths from different modes of burning may be 
difficult to distinguish. Therefore, two inflow velocities were used at each Damkohler number, 
specifically, /3 = 1 and /3 = 2. This means that if the turbulent flame burns inductively, the 
turbulent flame widths will differ by a factor of approximately two. However, if the flames 
are driven by turbulent mixing, then the turbulent flame widths should be independent of /3. 
This means that the crucial comparison required to determine the burning mode is between 
flames at different inflow speeds at the same Damkohler number, and direct measurements 
of local turbulent flame widths do not need to be evaluated, nor are comparisons of these 
widths necessary at different Damkohler numbers. 

Simulations of one-dimensional zero-turbulence inductive flames and three-dimensional 
inert turbulence were first obtained, and then superimposed to initialize each calculation. 
Each simulation was run with a resolution of 256 x 256 x 1024. Adaptive mesh refinement 
was not used. Table [[] gives the conditions for the six simulations. 



4. RESULTS 

Figures [3]^a,b) show two-dimensional vertical slices (through three-dimensional simula- 
tions) of burning rate (left) and temperature (right) for the Dai6 = 1/3 cases, (a) (3 = 1 and 
(b) (5 = 2. The two central panels in each figure show a snapshot of the turbulent simula- 
tions, and the narrow edge panels show the corresponding images for the (zero turbulence) 
inductive flames for comparison. The white lines show three relevant length scales (/, Zi6, 
and Ai6). Note that only Zi6 differs between the two cases (due to the dependence on u$ and 
therefore (3). Warm oxygen is being fed from below at uq = (3u, burns to sulphur, and leaves 
the domain through the top boundary. In both cases, there is a large volume of fuel burning, 
many times the integral length scale. It is clear that there is a high level of turbulent mixing, 
but the width of each flame appears to be roughly the same as the laminar inductive flame at 
the corresponding inflow speed. Importantly, the (3 = 2 flame appears significantly broader 
than for (3 = 1. This suggests that for Dai6 ~1 5 the flames burn inductively. 

Figures E£a,b) show the corresponding slices for the Dai 6 = 3 cases (the Dai 6 = 1 
cases present intermediate behavior, and are not shown). Note the domain size and integral 
length scale are 27 times larger than figure [3j and the inflow and turbulent intensity is 3 times 
higher. Correspondingly, the inductive length Iiq and A-width are relatively shorter here. In 
both cases, the local turbulent flame width appears to be narrower than the corresponding 
inductive flames shown, and is even filamentary in places. The crucial point to note here 
is that the local turbulent flame width does not differ significantly between the (3 = 1 and 
(3 = 1 cases. Therefore, turbulent mixing must be driving the flame propagation. If the 
flames were burning inductively, the (3 = 1 case should be broader than the (3 = 1 case. This 
is evidence that oxygen burns as a A-flame for Dai6 > 1- 

An underlying assumption of the scaling analysis in section [2] is that the nuclear time 
scale Tie is constant. This assumption, combined with a turbulent flame width narrower 
than / 16 , suggests that the turbulent flame speed is faster than si 6 (i.e. u = (3u). Figure [5] 
shows the flame displacement £(£) = z (t) — z (0) from the initial position z (0) as a function 
of time for all six cases, where the flame position has been defined as 



*>(*) = , / y n / pXudV. (7) 

A {pX 16 ) J v 

The displacement has been normalized by the integral length scale /, and the time has been 
normalized by the integral length eddy turnover time r = l/u. It is clearly evident that 
the Da30Bl flame is indeed burning significantly faster than uq and propagates towards the 
inflow boundary; we note that it is the high Damkohler case (Da30Bl) that is of particular 
interest. 
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This is an interesting consequence of turbulence- driven oxygen flame propagation, be- 
cause it suggests that the conventional idea that the oxygen flame is locked at some distance 



(depending solely o n uq and r 16 ) behind the carbon flame (e.g. iTimmes fc Woosleyi 11992 



Lisewski et al.ll2000l ) cannot be the case under these conditions. For an order-of- magnitude 
analysis, assume the turbulent burning time scale of carbon at p = 2.5 x 10 7 is of the or- 
der of 10~ 3 s, which means the carbon A-flame local speed and width are approximately 
s^ 2 rj 3 x 10 6 cm/s and A12 ~ 3 x 10 3 cm, respectively. The corresponding values for the 
resulting oxygen flame are s^ 6 = 1.13 x 10 7 cm/s and Ai6 = 1.81 x 10 5 cm, respectively. The 
oxygen A-flame is about an order of magnitude faster and two orders of magnitude thicker 
than the carbon flame. This suggests that the flame actually burns as a single compound 
carbon-oxygen flame with local flame speed and width close to that of the oxygen A-flame. 

Figure M shows two-dimensional slices through a three-dimensional simulation of a com- 
pound carbon-oxygen flame where the inflow velocity and turbulent intensity were matched 
to case Da30Bl. We emphasize that at this resolution, the carbon flame is far from being 
well-resolved. The simulation was initialized with a discontinuity halfway up the domain, 
with cold fuel under hot ash. Specifically, the fuel was at a density and temperature of 
p = 2.5 x 10 7 g/cm 3 and T = 6 x 10 8 K, and consisted of 40% carbon and 60% oxygen. The 
initial ash consisted of 40% magnesium and 60% sulphur, with density and temperature of 
p fa 1.39 x 10 7 g/cm 3 and T « 3.8 x 10 9 K, respectively, and was allowed to evolve to the 
appropriate state. The figure panels are carbon mass fraction, oxygen mass fraction, oxygen 
burning rate and temperature, respectively. It appears that, as expected, turbulent mixing 
is able to drive a compound flame. 



The potential transition to detonation suggested by IWoosley et al.l (120 lOl ) requires the 
formation of a region of approximately 10 km at a temperature of approximately 3.6 x 10 9 K. 
To investigate the existence of such a region, the temperature field from the DalOB2 case 
averaged using a top-hat cubic filter of size 33 3 cells, which corresponds to a length for each 
side of approximately 2.33 km. The filtered temperature and temperature variance were 
found at each point in space and time, and are plotted in the form of a joint probability 
density function in figure [7J The red curve denotes the zero turbulence case filtered in the 
same way. The yellow line denotes the minimum variance achieved for each temperature in 
the range denote d by vertical black lin es, over which a potential transition to detonation 



was proposed by IWoosley et al.l (120101 ) . Almost every point within the flame lies above 
the laminar profile, but there are low probability events that have a low variance within 
the required temperature range. However, the minimum variance over this range (denoted 
by the green c i rcle) is approximately 3.5 x 10 15 K 2 . The candidate case from figure 3b of 



Woosley et al.l (120101 ) has a temperature of approximately 3.67 x 10 9 ±2.5 x 10 7 K. Assuming 
a uniform distribution gives a variance of approximately 2 x 10 14 K 2 , which is an order of 
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magnitude less than in the current simulation. Allowing for a larger range, say ±le8K, 
gives a variance of approximately 3 x 10 15 K 2 . These estimates are not much lower than 
the laminar flame (shown by the red curve), which suggests that the turbulence does not 
lead to the formation of a plateau or ledge under these conditions. Furthermore, the filter 
size here is much smaller than required by about a factor of 4, and the variance will only 
increase with a larger filter size. This does not provide support for the suggests transition to 
de tonation in oxygen fo r the conditions studied, but it should be noted that the simulations 



m 



Woosley et al.l ( 120101 ) had a turbulent power over 15 times that used here. 



5. CONCLUSIONS 

The theo r etical treatment of distributed carbon-burning thermonuclear flames from 



Asp den et al.l ( 120101 ) has been applied to the trailing oxygen flames and compared with 
three-dimensional simulations over a range of Damkohler numbers. It was shown that for 
Dai6^1, turbulence does not greatly alter the flame from one in which the oxygen burns 
purely inductively. Since turbulence accelerates the carbon flame however, the width of the 
oxygen flame is enormously broader than in the laminar case. For Da^^l, turbulence en- 
hances heat transfer and drives flame propagations that is narrower than the corresponding 
zero turbulence inductive oxygen flame. This is somewhat counterintuitive as turbulence 
typically broadens interfaces rather than sharpening them. A consequence of burning in 
this limit is that the oxygen can burn faster than the inductive flame speed (but is limited 
by the carbon flame speed of course). Therefore, the oxygen flame does not trail behind 
the carbon flame (at a distance equal to the post carbon-flame velocity times the oxygen 
burning time scale), but burns as a compound carbon-oxygen. This suggests that a single 
level set is a suitable flame model for the compound flame under these conditions. Averaging 
the temperature field using a cubic filter suggested that the temperature variance in at the 
desired con ditions is too h i gh to support the potential transition to detonation in oxygen 



proposed in IWoosley et al.l ( 120101 ). However, this does not preclude this kind of transition 
under different conditions, such as higher turbulence or lower densities, and it should be 
borne in mind that only one such event would be required, and in the star there are many 
realizations. 
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Fig. 1. — Scaling relations for oxygen flame widths as a function of Dai6- The red lines 
show the inductive flame width Z 16 (solid for (3 = 1, and dashed for (3 = 2). The blue line 
shows the turbulent width t[ 6 if turbulent mixing drives the flame, which scales with Dai6 
for Dai6 ^1 an d is equal to the A-flame thickness for Dai6 ^1 • 
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Fig. 2. — Diagram of the simulation setup (shown in two- dimensions for clarity). The domain 
is initialized with a turbulent flow and a flame is introduced into the domain, oriented to 
that the flame propagates downwards against the imposed mean flow. The turbulence is 
maintained by adding a forcing term to the momentum equations. The top and bottom 
boundaries are outflow and inflow, respectively. The side boundaries are periodic. 
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Case 


Da03Bl 


Da03B2 


DalOBl 


DalOB2 


Da30Bl 


Da30B2 


Damkohler number (Da^) 


1/3 


1/3 


1 


1 


3 


3 


Inflow factor (/3) 


1 


2 


1 


2 


1 


2 


Domain width (L) [km] 


3.50 


3.50 


18.1 


18.1 


100 


100 


Domain height (H) [km] 


14.0 


14.0 


72.4 


72.4 


400 


400 


Integral length scale (I) [km] 


0.350 


0.350 


1.81 


1.81 


10.0 


10.0 


Turbulent intensity (u) [km/s] 


65.4 


65.4 


113 


113 


200 


200 


Inflow velocity (u ) [km/s] 


61.46 


12.65 


10.05 


20.68 


15.73 


34.90 


Inflow density (p ) [xl0 7 g/cm 3 ] 


1.786 


1.739 


1.895 


1.846 


1.998 


1.940 


Inflow temperature (T ) [xl0 9 K] 


2.950 


3.035 


2.790 


2.870 


2.656 


2.740 



Table 1: Simulation properties. The f denotes the synthetic inflow conditions that were 
chosen to position the laminar inductive flame close to the half-way point in the domain. 
The desired conditions were po = 1-69 x 10 7 g/cm 3 , T = 3.14 x 10 9 K and Uq = f3u. 
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Fig. 3. — Two-dimensional vertical slices (through three-dimensional simulations) of burning 
rate (left) and temperature (right) for the Dai6 = 1/3 cases, (a) /3 = 1 and (b) (3 = 2. The 
two central panels show a snapshot of the turbulent simulations, and the narrow edge panels 
show the corresponding images for the (zero turbulence) inductive flames for comparison. 
The white lines show three relevant length scales (/, 1±q, and A^). 
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Fig. 4. — As figure [3] for the Dai 6 = 3 cases. 



-18- 






03 A 
o -4 



CL 
03 



-5 



-9 



-10 







Da03Beta1 
Da03Beta2 
Da10Beta1 
Da10Beta2 
Da30Beta1 
Da30Beta2 



10 



15 



20 



25 



30 



35 



Time 



Fig. 5. — Flame displacement from the initial position as a function of time for all six cases. 
The displacement has been normalized by the integral length scale /, and the time has been 
normalized by the integral length eddy turnover time r = l/u. 
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Fig. 6. — Two-dimensional slices through a three-dimensional simulation of a compound 
carbon-oxygen flame. The panels are carbon mass fraction, oxygen mass fraction, oxygen 
burning rate and temperature, respectively. 
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Fig. 7. — Joint probability density function of temperature and temperature variance evalu- 
ated using a 2.33 km cubic filter for case DalOB2. The red curve denotes the zero turbulence 
case filtered in the same way. The yellow line denotes the minimum variance achieved for 
each temperature in the range denote d by vertical bl a ck lin es, over which a potential tran- 
sition to detonation was proposed by IWoosley et al.l (|2010|). The green circle denotes the 
lowest variation that was found within this temperature range. The inset shows a zoom of 
the conditions of interest. 



